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We study the theory of itinerant-hole photoluminescence of two-dimensional electron systems in 
the regime of the magnetically induced Wigner crystal. We show that the exciton recombination 
transition develops structure related to the presence of the Wigner crystal. The form of this structure 
depends strongly on the separation d between the photo-excited hole and the plane of the two- 
dimensional electron gas. When d is small compared to the magnetic length, additional peaks 
appear in the spectrum due to the recombination of exciton states with wavevectors equal to the 
reciprocal lattice vectors of the crystal. For d larger than the magnetic length, the exciton becomes 
strongly confined to an interstitial site of the lattice, and the structure in the spectrum reflects 
the short-range correlations of the Wigner crystal. We derive expressions for the energies and the 
radiative lifetimes of the states contributing to photoluminescence, and discuss how the results of 
our analysis compare with experimental observations. 

PACS numbers; 73.20.Dx, 71.35. -l-z, 78.20.Ls 
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I. INTRODUCTION 

Recent experimental studies of high-mobility 
GaAs/AlGaAs devices using itinerant-hole photolumi- 
nescence have shown this to be a sensitive probe of the 
two-dimensional electron systems formed in these devices 
. Features in the photoluminescence spectra have been 
found to correlate well with transport measurements of 
the integer and fractional quantum Hall effects 
and the insulating regime observed at very low filling 
fraction [^-0. These studies are particularly interesting 
in this low filling fraction regime, where transport mea- 
surements can provide only limited information on the 
state of the two-dimensional electron gas. In particular, 
they may provide information on the transition to the 
Wigner crystal which is the expected groundstate under 
these conditions. Related experimental studies of the 
photoluminescence of intentionally acceptor-doped sam- 
ples, in which the photo-excited hole becomes bound to 
an acceptor, have also found structure associated with 
the transition to an insulating regime |^,^. However, the 
signature of the insulating state in those experiments is 
quite different from that in the undoped devices in which 
the photo-excited hole is free to move ("itinerant-hole" 
photoluminescence) . 

In the undoped samples, the transition to the insulat- 
ing regime is associated with the appearance of additional 
lines in the photoluminescence spectrum (5|-|^]. It has 
been suggested that the appearance of these lines does 
indeed indicate the formation of a magnetically induced 
Wigner crystal |^,||. However, the theory underlying 
itinerant-hole photoluminescence in the extreme quan- 
tum limit is not well understood. Existing theories of 
photoluminescence in the Wigner crystal regime apply 
only to the limit in which the photo-excited hole is far 
from the electron gas compared to the typical electron- 
electron spacing pn,pd|. This condition is not appropri- 



ate for typical devices in which itinerant-hole photolu- 
minescence is observed: while the asymmetry of these 
devices does cause the hole to lie some distance from the 
electron gas, it is believed to be close compared to the 
electron-electron spacing . 

We present a theory for itinerant-hole photolumines- 
cence in the Wigner crystal regime which applies to sys- 
tems in which the photo-excited hole is close to the two- 
dimensional electron gas. We study the limit of strong 
magnetic field in which all the electrons and the photo- 
excited hole are confined to states in the lowest Lan- 
dau level. To represent the asymmetry of the quantum 
well, the hole is assumed to lie in a plane a distance 
d away from the two-dimensional electron gas (which is 
also assumed to have zero thickness). Such a model is 
common in theories of photoluminescence in the frac- 
tional quantum Hall regime |]l3|-|l5|]. In that regime, the 
model presents an intrinsically strongly-coupled many- 
body problem, which, for the most part, has required 
numerical investigation. For very low filling fraction and 
for d small compared to the typical electron-electron 
spacing, certain simplifications arise. Under these con- 
ditions, we expect that the low-lying energy states in the 
presence of the photo-excited hole involve the formation 
of an exciton, with one electron strongly bound to the 
photo-excited hole. A well-defined exciton can form pro- 
vided the size of the exciton, set by the magnetic length 
I = y^h/eB, is small compared to the electron-electron 
spacing. Moreover, since the exciton is neutral, it will 
couple rather weakly with the other electrons. In particu- 
lar, we expect that for small d, the presence of the exciton 
will not significantly perturb the groundstate (Wigner 
crystal) configuration of the other electrons, and the ex- 
citon will behave as a rather non-invasive probe of this 
crystalline state. This expectation is motivated by stud- 
ies of the classical groundstate of a system of electrons in 
the presence of an ionized donor impurity located a dis- 
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tance d from the plane of the two-dimensional electron 
gas It is found that for d < 0.29a (where a is the lat- 
tice constant of the Wigner crystal) the system adopts a 
groundstate configuration in which one electron is bound 
to the donor impurity, while the remaining electrons lie 
close to the sites of a triangular lattice. 

In this paper, we discuss how the photoluminescence 
spectrum arising from the recombination of the exciton 
is affected by the presence of the crystal. We study a 
model in which the exciton moves in the static potential 
set up by a triangular lattice of electrons. Our analysis 
neglects the dynamical properties of the lattice (phonon 
coupling). This approximation is justified, as the struc- 
ture that we will find is on a energy scale large compared 
to the typical magnetophonon energy, 0{Pe^ /ea^): the 
lattice can therefore be considered to be static on the 
timescale necessary to define these exciton states. In 
the present work we neglect exchange processes between 
the electron in the exciton and those forming the lattice; 
these will be considered in a separate paper [0 . 

The behaviour of our model is strongly dependent on 
the separation d. For small d, the coupling of the exci- 
ton to the lattice is weak; in this case, additional peaks 
appear to higher energy of the main exciton line arising 
from the recombination of exciton states with momenta 
equal to reciprocal lattice vectors of the crystal. The 
spectrum effectively carries diffraction information on the 
crystal and is therefore sensitive to long-range crystalline 
order. As d increases, the coupling becomes stronger and 
the low-lying exciton states become increasingly confined 
to the interstitial regions of the crystal. These states are 
less sensitive to the presence of crystalline order, but de- 
pend on the short-range correlations of the electron gas. 
The crossover from weak to strong coupling occurs when 
the (dipole) potential energy of the exciton, 0{d^e^ /ea'^), 
becomes larger than the kinetic energy cost to confine the 
exciton to a size of order the lattice constant /Ma^ 
(where M is the effective mass of the exciton). We show 
that this crossover occurs when the separation between 
electron and hole planes d becomes larger than the mag- 
netic length. 

We compare the results of our model with experimental 
observations of photoluminescence in the extreme quan- 
tum limit. We find that the energy scale predicted by our 
model for the splitting of the exciton transition is com- 
parable to the energy scale of the structure observed in 
experiment. We argue that the observed peak-splitting 
that has been associated with the Wigner crystal regime 
could arise from the groundstate and first excited state 
of a strongly-confined interstitial exciton. However, other 
explanations for this splitting cannot be ruled out, and 



it may be that the structure that we predict is not seen 
due to a lack of population of the higher energy exciton 
states. Further experimental work is required in order to 
identify this structure. We argue that this could best be 
observed in optical absorption experiments. 

The paper is organized as follows. In the following 
section we review the theory of the free two-dimensional 
exciton in a strong magnetic field. We derive expressions 
for the binding energy and effective mass of the exciton 
for the case in which the electron and hole lie in planes 



separated by d. In section [II, we study the coupling of 



the exciton to the periodic potential set up by a Wigner 
crystal of other electrons. We derive exact expressions 
for the matrix elements of the periodic potential of the 
crystalline lattice within the basis of free exciton states, 
and study the resulting low-energy exciton states in the 
weak-coupling limit, d <C ^, using perturbation theory 
about the free exciton states. In section [V we study the 
limit of strong coupling, d ^ I. To do so, we develop an 
effective Hamiltonian for the motion of the exciton in a 
smooth external potential. We apply this to the case in 
which the potential is due to the presence of a Wigner 
crystal, and study the low-energy eigenstates close to an 
interstitial site of the electron crystal. For both the weak 
and strong coupling limits we derive expressions for the 
energies and relative radiative lifetimes of the low-lying 
energy states. In section ^ we discuss how these results 
compare with existing photoluminescence measurements 
on GaAs/AlGaAs heterostructures in the extreme quan- 
tum limit. Finally, section VI contains a summary of the 
main points of the paper. 



II. FREE TWO-DIMENSIONAL EXCITON 
STATES IN STRONG MAGNETIC FIELD 

The wavefunctions and the energy spectrum of a free 
two-dimensional exciton in a strong magnetic field have 
been discussed by Lerner and Lozovik |18|]. T hey built on 
the work of Gor'kov and Dzyaloshinskii |19[| who showed 
that, within the effective mass approximation for the con- 
duction and valence bands, the exciton states can be de- 
scribed in terms of a conserved two-dimensional momen- 
tum, P. In the limit of strong magnetic field, when the 
cyclotron energies of both the electron and hole become 
large compared to the typical electron-hole interaction 
energy, the exciton wavefunction is completely specified 
by the momentum P and the Landau level indices of the 
electron and hole. The lowest energy exciton state, in 
which the electron and hole are both in the lowest Lan- 
dau level, has the wavefunction [n9[ 
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where and are the electron and hole positions, 
rp = z X PP/h, A is the area of the system, and 
the vector potential has been chosen in the symmetric 
gauge, A{r) = B x r/2. We have chosen units in which 
h = I = 1, and, in the following, we express energies in 
units of e^/Aneeol- To make the discussion more trans- 
parent, however, we reintroduce these units at appropri- 
ate points. 

The energies of the free exciton states, relative to the 
zero-point kinetic energy of the electron and hole, are 
given by the expectation values of the electron-hole in- 
teraction potential We study a situation in which 
the electron and hole move in planes separated by a dis- 
tance d, such that the interaction is 



1 



47reeo 



(2) 



The interaction energy of the exciton state \P) is there- 
fore 
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where Jq{z) is the ordinary Bessel function. We have 
not found a closed-form expression for this integral (in 
the limit d = it reduces to the expression derived in 
Ref. However, in what follows, we are interested in 

exciton states with wavevectors on the scale of the recip- 
rocal lattice vectors of the crystal, 0(1/ a). At low filling 
fraction, this is much smaller than the scale on which 
the energy (^) varies, 0{\/l), so for our purposes it is 
sufficient to work with an expansion 
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(4) 



which may be interpreted in terms of a "binding energy" , 
and an "effective mass", Af^, 
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.J^- d+l/2.J^(f + 0{cP), d<^l (5a) 

1/^ + 0(1/^2), d> 1 (5b) 

e'^^'^{l + d^) erfc [d/Vfj - d/2 (6) 

d + ?,/2^/^d^ + 0{d^), d-^l (6a) 
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In many ways the exciton behaves as a simple free par- 
ticle with a mass M^. In particular, the velocity of an 



exciton in the state \P) is dE{P)/dP In view of 

this, it is convenient to think of the contribution to Ed{P) 
which is momentum-dependent as the "kinetic energy" of 
the exciton, even though this energy originates from the 
electron-hole interaction. The exciton is overall charge 
neutral, but it does carry a dipole moment and therefore 
couples weakly to an external electrostatic potential. The 
dipole moment is of size —edz perpendicular to the plane 
of the interface, and of size —erp parallel to this plane. 
It is the dipole moment of the exciton that will cause 
it to be scattered by the Wigner crystal of other elec- 
trons. The strength of the coupling of the exciton to a 
Wigner crystal with a lattice constant a is determined by 
the competition between the typical dipole energy of the 
exciton, d^e^/ea'^, which is minimized at an interstitial 
site of the lattice, and the kinetic energy cost to confine 
the exciton to such a site, l/{A'Ida'^)- Using the asymp- 
totic expression for the mass (6b), it is found that these 
two energies become equal when d/l ~ l/v^^^'^, where 
h' is the filling fraction of the two-dimensional electron 
gas. This condition is so weakly dependent on the filling 
fraction that it is accurate to say that one expects the 
exciton states to crossover from being weakly coupled to 
the lattice when d < Z to strongly confined to interstitial 
regions when d^ I. 

In order to understand the photoluminescence spec- 
trum it is essential to know both the energies of the ex- 
citon states and the rate at which they decay to emit 
radiation. Within the effective mass approximation for 
the electron and hole bands, the radiative lifetime of 
an exciton state depends on two factors: the band-to- 
band dipole matrix element of the host semiconductor, 
and a matrix element between the electron and hole en- 
velope functions. The first contribution is constant for 
all transitions. We study only the envelope term, which 
is sufficient to describe the relative recombination rates 
of different exciton states. The operator that describes 
electron-hole annihilation may be written in second quan- 
tized notation as 



L = 
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d'^r ■(/ie(r)i/)ft(r). 
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where iph and ipe are the field annihilation operators for 
holes and electrons, and the wavevector of the emitted 
photon has been assumed to be zero. We have chosen 
a convenient normalization, which is arbitrary since we 
consider only relative radiative decay rates. We define 
the "luminescence strength" of a transition between any 
pair of (many-body) states |i) \f) by 



L.f = \{f\m\' 



(8) 
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For the free exciton states discussed above (|l|) the only 
available final state is the vacuum and we find 



Lp = \{YAC\L\P)\^ =6p,„. 



(9) 



Thus, due to momentum conservation, only the P = 
exciton state can emit a long-wavelength photon. Since 
the free exciton states form a complete set of states for 
the exciton in the lowest Landau level, the luminescence 
strength of a transition in which an exciton in a general 
state \ip) decays to emit radiation may be found from 
Eqn. (|9|) to be given by 



(10) 



We will use this expression in the following discussions 
which concern the exciton states in an external potential, 
for which the momentum is not a good quantum number. 

The relative intensity of a photoluminescence transi- 
tion is proportional to its luminescence strength multi- 
plied by the probability for a photo-excited hole to be 
in the initial state of the transition. These probabili- 
ties are difficult to quantify: the lifetime of the photo- 
excited hole can be shorter than its equilibration time, 
so one may observe recombination from high-energy non- 
equilibrium states. Although we are primarily interested 
in photoluminescence, we note here that the lumines- 
cence strength Li f also characterizes the strength of the 
transition, as observed in optical absorption, from the 
state I/) to the state \i). In this case, uncertainties re- 
lated to non-equilibrium populations do not arise as the 
lifetime of the hole does not limit the time available for 
the initial states to equilibrate. 



III. EXCITON STATES IN THE PRESENCE OF A 
WIGNER CRYSTAL: WEAK COUPLING LIMIT 



We will now study the scattering of the exciton 
states (|l|) by the electrons forming the Wigner crystal. 
As described in the introduction, we neglect the exchange 
interaction between the electron in the exciton and those 
forming the crystal. We anticipate that this is a good ap- 
proximation at low filling fraction, when the fraction of 
basis states excluded from the exciton wavefunction by 
the Pauli exclusion principle is small; this is confirmed 
by numerical studies in which the effects of exchange are 
included |l^. The model we discuss also correctly rep- 
resents a situation in which the spin or subband index 
of the electron in the exciton is different from the cor- 
responding index of the electrons in the two-dimensional 
electron gas. 

Neglecting exchange interactions, the exciton is scat- 
tered only by the charge density 



2nP 



which represents a system of electrons in lowest Landau 
level orbitals centered at the sites {Ra} of a triangu- 
lar lattice with a lattice constant a (the filling fraction 
is therefore ly — {An / ^/3)P / a^) . The magnitude of p is 
chosen to provide a uniform neutralizing positive back- 
ground. The resulting external potential energies of the 
electron and hole are 
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The motion of the exciton in this potential is fully described by the matrix elements of all interactions within the 
basis of free exciton states (^ 



(13) 



The first term is diagonal in this basis and gives rise to the free exciton energy (y,^) ■ The last two terms describe the 
scattering. Due to the discrete translational symmetry of the crystal, the only matrix elements which are non-zero 
are those between states which differ by a reciprocal lattice vector of the crystal, {K}. Since we are ultimately 
interested in states that can emit radiation, and which therefore from Eqn. (10) must contain some component of the 
zero-momentum free exciton state, we need only study the basis states \K) in which the wavevector of the exciton is 
a reciprocal lattice vector. After a lengthy but straightforward process wc find 



Hk'.k = Sk'.k Ed{K) + (1 - 5K',K)y- 



-■i(K-K'f /i 

\K-K'\ 



^iKy.K'.z/2 _ ^-\K-K'\d^-iKxK' .z/2 



(14) 



While the above expressions for the matrix elements are exact and therefore fully describe the behaviour of our 
model, in their present form they are not useful for analytic purposes. To obtain some insight into the properties of 
this system we therefore introduce some approximations. In what follows, we develop a perturbation theory in d/l, 
based on the free exciton states, which is valid in the limit of low filling fraction. 
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Consider the case d = in which the hole moves in 
the same plane as the electrons. In this limit, the zero- 
momentum exciton state is not scattered by the poten- 
tial, and is therefore an exact energy eigenstate 

Hk',k=o = Ed=o{Q) 6k',q. (15) 

All other energy eigenstates must be orthogonal to this 
state, so from Eqn. ( p^ they must have zero lumines- 
cence strength. There is therefore only a single line in 
the photoluminescence spectrum. Moreover, this line 
appears at the same energy as the exciton line in the 
absence of the two-dimensional electron gas. This recov- 
ers a restricted version of the general result that when 
electrons and holes are confined to the lowest Landau 
level and move in the same plane, photoluminescence 
contains no spectroscopic information on the state of the 
two-dimensional electron gas ||2^,|l^ . 



Although the free exciton state with zero momentum is 
an exact energy eigenstate of the system at d = 0, other 
free exciton states are not. There do remain non-zero off- 
diagonal matrix elements coupling these states. However, 
in the limit of small filling fraction, this coupling may be 
neglected for the lowest energy states: states with small 
momenta \K\ ^ 1/a, are only coupled by off-diagonal 
terms of order (l/a)^, whereas the energy spacing be- 
tween these states is of order fi^ /{Md=oo?) ~ (l/a)'^- We 
choose to work in an expansion in small filling fraction, 
and neglect terms of order {l/aY and higher. Within this 
approximation, the free exciton states are eigenstates of 
the Hamiltonian at d = 0. 

We now study the deviations from d = within a per- 
turbation expansion in the separation d. Expanding the 
Hamiltonian in this parameter, we find 



J 



Hk' 



K 



TK'K 



+ (1 - 5K'M)vd + 0(l/a3) + O(dVa')- 



(16) 



To lowest order in d all states are coupled by a matrix 
element (1 — 5k' .K)i^d. This leads to a mixing of the ra- 
diative state, \K = 0), with all other exciton states of re- 
ciprocal lattice vector; if these states become populated 
on photo-excitation, each will contribute an additional 
line to the photoluminescence spectrum. The strongest 
new line arises from the six lowest-lying energy states, 
at \K\ = 47r/(V3a). Within degenerate perturbation 
theory, these states are split into a five-fold degenerate 
level shifted in energy by —vd, and a single state shifted 
by -\-bi'd. We now apply first-order perturbation the- 
ory to calculate the mixing oi\K = Q) into these states. 
This amplitude is zero for the five degenerate states, and 
equal to ^/^ud for the single state split off from these five. 
The luminescence spectrum therefore consists of a strong 
spectral line from what was the free zero-momentum exci- 
ton state, at an energy (relative to the band-gap plus zero 
point kinetic energy) and with a luminescence strength 



Ea — —Bd, 

io = 1 - O(rf') 



(17) 
(18) 



and an additional line arising from the exciton states with 
the smallest non-zero reciprocal lattice vector, with an 
energy and luminescence strength 



El 



u = 



-Bd 



(47r/\/3a)2 



f(4^/V3a)2 



^'4d^ 



bvd (19) 
(20) 



The next line to higher energy arises from the six states 
at \K\ = An /a, and has a luminescence strength that is 
smaller than this by a factor of [(47r/V3)/(47r)]'^ = 1/9. 



Thus, for small d the presence of the crystal causes 
peaks to appear in the photoluminescence spectrum to 
higher energy than the main exciton line. The energy 
spacing between these peaks is primarily determined by 
the effective mass of the exciton and the reciprocal lat- 
tice vectors of the Wigner crystal, with a small correction 
due to the coupling to the lattice potential. The lumines- 
cence strengths of these higher energy exciton states are 
small compared to the lowest energy exciton transition, 
and will therefore give rise to much weaker spectral fea- 
tures. The luminescence strengths of these features grow 
as d increases. However, at the same time, the exciton 
becomes more strongly coupled to the lattice potential, 
so the energies of these states deviate from the energies 
of the reciprocal lattice vector excitons. Once the cou- 
pling is larger than the spacing between the free exciton 
states, i/d ^ l/M^a^, the corrections to these energies 
become large, and the above perturbation theory breaks 
down. This condition is equivalent to the strong-coupling 
condition, d ^ I, that we derived earlier. In the strong- 
coupling limit, a quite different approach is required to 
describe the low energy exciton states. 



IV. EXCITON STATES IN THE STRONG 
COUPLING LIMIT: THE INTERSTITIAL 
EXCITON 



We now turn to the limit of strong exciton-lattice cou- 
pling, in which the kinetic energy of the exciton is small 
compared to its dipole potential energy in the lattice. If 
we were to neglect the kinetic energy completely, then 
the groundstate of the system would be the same as that 
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of the classical groundstate of electrons in the presence 
of an ionized donor impurity the exciton, with its 
weak dipole moment, would position itself at an intersti- 
tial site of the lattice, and, for d < 0.29a, the remaining 
electrons would form an essentially undeformed Wigner 
crystal. The exciton does, however, have some residual 
dynamics. In this section we develop an effective Hamil- 



tonian that describes the motion of the centre of the ex- 
citon in a smooth external potential. We then apply this 
to the potential close to an interstitial site of the Wigner 
crystal (|l2|), and study the energies and luminescence 
strengths of the low energy states. Those who do not 
care to follow the derivation of the effective theory may 



proceed to section IV B without loss of continuity. 



A. Effective Hamiltonian 



In order to derive an effective Hamiltonian for the motion of an exciton in a smooth external potential, we extend 
the approach of Gor'kov and Dzyaloshinskii to include the potential energies V'^{re) and V'^{rfi) felt by the 
electron and hole. We start from the full Hamiltonian for the interacting electron-hole pair in a uniform magnetic 
field 

and work in the symmetric gauge. We transform to a new set of co-ordinates, the first of which is the momentum 
defined in Ref . [ p^ ; the remaining co-ordinates are the position of the centre of mass and a momentum and position 
describing the internal motion 

P = Pe+Ph-^^Bx{re-rh) (22) 

Rc = — (23) 

rUe + urih 

rrie + mh 2 \ nie + nih J 
r = re-rh. (25) 

The centre-of-mass and the internal co-ordinates behave as independent, canonically conjugate pairs. In terms of 
these new co-ordinates, the Hamiltonian is 

^ = ^^^H^ + iP + iel2B.rf ^ ^ ^ ^^^^^ ^ _ ^^^^^ ^^^^ 



where M = mg + m/i, /i = merrihl {nie + ?7i/i), 7 = 
{mh - me)/ime + mu), and rja = ma/ (me + ruh). 

First consider the free exciton, — — 0. The 
Hamiltonian is independent of Rc, so the momentum, 
P, is conserved and can be replaced by its eigenvalue to 
leave a Hamiltonian for the internal co-ordinates, r and 
p. If the interaction potential, V'^'^ {r), is neglected in re- 
lation to the kinetic energy, the internal motion consists 
of Landau level orbits centred on rp = z x PP /Ti. This 
is the neglect of Landau level coupling, and is the simpli- 
fication used in Refs. |p^ , |l9| ; in the lowest Landau level, 
it leads to the states that we have discussed above (|l|). 

We now introduce the potentials, and V'^, and fol- 
low a similar approximation. The momentum is no longer 
conserved. However, provided the potential is sufhcicntly 
weak, the momentum only changes slowly compared to 
the rapid cyclotron motion. The motion can then be 



treated adiabatically, with the fast internal motion fixed 
in the lowest Landau level and adjusting to follow the 
slowly changing momentum. This procedure is analogous 
to the Born-Oppenheimer approximation in the theory 
of molecular dynamics |2^ , where the fast electronic de- 
grees of freedom are eliminated to provide an effective 
theory for the slow atomic co-ordinates. In the same 
way, we obtain an effective Hamiltonian for the opera- 
tors P and Rc- Before we write down this Hamiltonian, 
we make one last change of variables. It is convenient to 
work is terms of the co-ordinate R = {r^ + rh)/2 rather 
than the centre of mass. This avoids any irrelevant mass 
dependence in our analysis ( "irrelevant" since both elec- 
tron and hole are restricted to the lowest Landau level). 
This transformation does not affect the commutation re- 
lation of the position and momentum co-ordinates. The 
effective Hamiltonian in terms of these co-ordinates is 
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H = + Ed{P) + j [V\R + (r' + rp)/2) + - (r' + rp)/2)] ^'r', 



(27) 



where hu>a/2 — heB/2fi is the zero-point kinetic energy of the electron and hole, and Ed{P) is the dispersion 
relation (0) arising from their mutual attraction. 

To convert this expression into a more convenient form, we make use of the simplifications available in the problem. 
As we are interested in small momenta, 0{l/a) <C l/l, we expand the dispersion relation to quadratic order. Also, 
provided the exciton is not very close to a lattice site, the potentials V and are smooth on the lengthscale of the 
magnetic length, and it is a good approximation to expand these terms in the integral of Eqn.(p7|) to first order in 
(r' + rp). We find 

H ~ hcj^/2 -Bd + ^^+ V^iR) + V\R) - ip. [VViR) - VV\R)] x z. (28) 

This is a general expression for the motion of the exciton in smooth external potentials, and in a strong magnetic 
field. An expression for the mass Md was calculated earlier (^. Since the position and momentum operators are 
canonically conjugate co-ordinates, the energy eigenstates of the exciton follow from the solutions of the Schrodinger 
equation with the Hamiltonian ( p8| ) ; it is to be understood that the final term of this expression is symmetrized in the 
position and momentum operators, such that the Hamiltonian is Hcrmitian. Two approximations were used to derive 
this expression. Firstly, the rapid cyclotron motion was assumed to adiabatically follow the changing momentum; this 
is valid provided the spacing between the energy levels arising from the centre-of-mass motion is small compared to 
the cyclotron energy, tkij^ (no Landau level coupling). Secondly, we expanded the external potential to first order in 
(r' + rp), neglecting terms of order PW^V. 



B. Motion in the Wigner Crystal: Harmonic Approximation 

We now use the formalism developed in the previous subsection to study the motion of the exciton in a Wigner 
crystal. We therefore introduce the electrostatic interactions ( p^ ) as the external potentials in the effective Hamilto- 
nian (p8|). This expression for the effective Hamiltonian requires the external potential to vary slowly on the scale of 
the exciton size, I. We therefore expect this formalism to provide an appropriate description of the low-energy states 
when d I, in which the exciton remains far from the interstitial sites. However, we still assume that d/a <ti 1, and 
expand the effective Hamiltonian in this parameter to find 

~2M'd^8^^\R-Rc,\^~^ATTeeQ h ' \R- R^\3' ^ ^ 



where the sums run over the triangular lattice sites 
{Ra}, and the final term should be interpreted as the 
symmetrized product of the momentum-dependent and 
position-dependent factors. We have dropped the con- 
stant binding energy, —Bd, and the zero-point kinetic 
energy, ?iWp/2, and have also neglected a constant en- 
ergy, —vd^ which arises from the dipole moment of the 
exciton in the electric field of the neutralizing positive 
background. Each of the terms appearing in this effective 
Hamiltonian has a simple intuitive interpretation. The 
first two terms simply represent the kinetic energy of the 
exciton and the dipole energy for the exciton to be cen- 
tred at the position R. The final term is less familiar, but 
also arises from a dipole energy: in this case, due to the 
in-plane dipole moment of the exciton, which is of a size 
—erp for an exciton with momentum P. The resulting 
contribution to the Hamiltonian resembles the first-order 
coupling of a charged particle to a vector potential. 



I 

We have argued that for large d/l the exciton is con- 
fined at the interstitial sites of the triangular lattice. We 
now study the low-lying states of this effective Hamil- 
tonian within the simplest approximation, in which the 
potential is expanded to harmonic order about an inter- 
stitial site. Keeping only the contributions to the con- 
finement potential arising from the nearest three lattice 
electrons, we find 

P^ 12J3V3d^ - ^4 (i^ X P) z, (30) 
2Md 2 4 a5 2 a3 ^ ^ ' V ; 

plus a constant energy of 11.57d^ /a'^ due to the dipole 
energy at the interstitial site of the triangular lattice . 
We have again set e^/Aireeol = h = I = 1. Due to the 
rotational invariance of this Hamiltonian, its eigenstates 
may be classified by their angular momenta, R x P.z. 
We restrict attention to states with zero angular momen- 
tum, since only these states can emit long-wavelength 
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radiation. The zero angular momentum states of the two- 
dimensional harmonic oscillator have the wavefunctions 
and energies 



En = hn{2n + 1) 



where 



=2 



/ 243V3 
V 4 ^ 47reeoZ ' 

/ 4 47reeo ?i^a^ 
243\/3 J 



1/4 



(31) 
(32) 

(33) 
(33a) 

(34) 
(34a) 



and Cn{z) are the Laguerre polynomials. The expres- 
sions for d ^ I follow from the asymptotic expansion of 
the effective mass (6b). Note that the two-dimensional 



harmonic oscillator states with zero angular momentum 
are spaced in energy by twice the oscillator quantum, hft. 

Radiative recombination occurs from each of these 
states. For consistency with the definition of the lumines- 
cence strength in the weak coupling limit, we construct 
states consisting of a superposition of harmonic oscilla- 
tor states, placed at each interstitial site of a crystal with 
an area A. The only such linear combination with non- 
zero luminescence strength is the state in which all states 
are combined with the same amplitude and phase. From 
Eqn.(H) we find that the luminescence strength of such a 
state constructed from tp^" is 



V3a2 



V3a2 ■ 



Note that this is independent of the index n, so the 
ground state and all excited states have the same lumi- 
nescence strengths Such behaviour contrasts with 
the weak-coupling limit for which the lowest energy exci- 
ton state is much more strongly coupled to radiation than 
the higher energy states. As the spacing d is increased 
from zero, so that the system evolves from the weak- 
coupling to the strong-coupling limit, the higher energy 
exciton transitions grow in luminescence strength from 
zero (optically inactive) to eventually attain the same 
strength as the lowest-energy exciton state. 

In summary, within the harmonic approximation for 
the potential at an interstitial site of the lattice, the exci- 
ton contribution to the photoluminescence spectrum con- 
sists of a series of lines which are uniformly spaced and 
have relative intensities given by the relative populations 



of the states. This approximation is valid provided the 
typical size of the exciton state is small compared to the 
lattice constant. From equation (jsj), one finds that this 
depends only very weakly on the filling fraction, and oc- 
curs when d^l. This is the same as the strong-coupling 
condition under which the perturbation expansion of the 
previous section failed. 



V. EXPERIMENTAL COMPARISONS 



Several groups have reported studies of the itinerant- 
hole photoluminescence spectra of the two-dimensional 
electron systems formed in high- mobility GaAs/AlGaAs 
devices [^^||^,^,^. Measurements on single heterojunc- 
tions and single quantum wells show qualitatively the 
same behaviour in photoluminescence: in the fractional 
quantum Hall regime two spectral lines are observed; as 
the filling fraction is reduced below about v — 1/6, the 
higher energy peak is found to split to form a doublet. 
The filling fraction below which this doublet structure 
appears correlates well with the filling fraction at which 
transport measurements show a transition to an insulat- 
ing state Q . This has motivated claims that the doublet 
is associated with the formation of a magnetically in- 
duced Wigner crystal 1^,^. However, the origin of this 
doublet is still not well understood. In this section, we 
discuss whether this doublet can be accounted for in 
terms of the structure that we predict. We will com- 
pare the energy splitting of the doublet with the energy 
difference between the two lowest-lying exciton states of 
our model that are optically active. 

For numerical comparisons of our theory with exper- 
iment, we will focus on the results of Goldys et al. 
presented in Ref. This paper reports studies 

of the photoluminescence spectrum of a high mobil- 
ity GaAs/AlGaAs heterostructure with density 3.2 x 

(35) 

10 cm . Specifically we will concentrate on a field of 
B = 12. 5T which is appropriate for Fig. 4(a) of Ref. 

for which the doublet structure, with a splitting of 
0.5meV, is well developed. Under these conditions, the 
filling fraction is i' = 0.1, corresponding to a magnetic 
length of I = 72A and a lattice constant of a = 6OOA for 
the triangular Wigner crystal. 

The most uncertain parameter that enters our model 
is the spacing d. This distance is expected to depend on 
the details of the band-bending in the vicinity of the in- 
terface, and may even vary with magnetic field jl^ . We 
will assume values oi d — 50A and d = lOOA as small 
and large estimates for this quantity. These are consis- 
tent with recent numerical studies of the binding of an 
exciton to the interface of a single-heterojunction p2[ . 
Note that these values are also consistent with the con- 
dition d < 0.29a which is required by our theory in order 
that the exciton does not strongly deform the crystal. 
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The smaller of these two separations, d = 50A, is 
less than the magnetic length, so we expect our weak- 
coupling formulae to be appropriate. In the weak cou- 
pling limit, the difference in energy between the main 
exciton line and the lowest-energy optically-active state 
is primarily determined by the free exciton energy at a 
wavevector equal to the smallest non-zero reciprocal lat- 
tice vector \K\ — 47r/(-\/3a). Taking the dielectric con- 
stant of GaAs to be e = 12.53, we find from equation (|^) 
that the energy spacing is l.SmeV. 

The larger spacing, d = lOOA, satisfies d > I, and we 
should therefore use our strong-coupling formulae, which 
state that the energy spacing between adjacent optically- 
active exciton states is 2hQ, where Ml is given by equa- 
tion (|33|). Using the appropriate value for the effective 
mass of the exciton , this evaluates to an energy spac- 
ing of O.TmeV dX d — lOOA. Under these conditions, 
the characteristic size of the groundstate wavefunction is 
i?o = 150A, which is smaller than, but comparable to, 
the distance to the closest saddle point of the crystal po- 
tential, a/-\/l2 = 170A. This indicates that, under these 
conditions, the harmonic approximation is at the limit of 
its validity. 

An exact calculation of the energy spacings predicted 
by our model for these two cases is likely to lead to 
energies that are slightly smaller than the above esti- 
mates deviations from the weak coupling limit will 
tend to push the energy spacing towards the (smaller) 
energy spacing one would obtain by using the strong- 
coupling theory; away from the strong-coupling limit, 
the use of the harmonic approximation will over-estimate 
the strength of the confinement potential and hence also 
overestimate the energy spacing. We also expect a small 
reduction of both energies due to the finite thicknesses 
of the electron and hole subband wavefunctions in a real 
device. Even allowing for these reductions, the resulting 
energy spacings remain larger than the expected zone- 
boundary magnetophonon energy, 0.2meV [2^. Our ne- 
glect of the dynamical properties of the Wigner crystal 
is therefore a consistent assumption. 

The energy separation between the lowest two opti- 
cally active states in our model is on the same scale as 
the observed energy splitting of the doublet, O.SmeV 
It is therefore possible that this structure is due to the 
splitting of the exciton peak as described by our model. 
Note that, if it is indeed correct to attribute this struc- 
ture to the predictions of our model, one must assume 
that the experimental system is in the strong-coupling 
limit, in which the exciton is strongly confined to an in- 
terstitial site. There are two important reasons for this, 
in addition to the fact that the experimental energy spac- 
ing (O.SmeV) compares more favourably with a spacing 



of d = lOOA than with smaller values. Firstly, it is ob- 
served that the intensities and the radiative lifetimes of 
the two peaks forming the doublet are similar The 
luminescence strengths of the two transitions must there- 
fore be comparable. In the weak coupling limit, all higher 
energy exciton states have a much smaller luminescence 
strength than the lowest energy exciton line; one must go 
to the strong-coupling limit before these become equal 
[see equation ( p5[ ) and the following discussion]. Sec- 
ondly, the observed energies of the two transitions form- 
ing the doublet are independent of the temperature of the 
substrate over the rang e 0.1 - 3K g. Since the Wigner 
crystal is expected to melt at a temperature of 0.3K at 
these densities |Q , one cannot attribute this structure to 
a nearly- free exciton state: the energy of this state would 
be sensitive to the changing long-range correlations at the 
melting transition. Instead, it is appropriate to attribute 
this structure to a strongly-confined interstitial exciton, 
which is only sensitive to the short-range correlations of 
the electron gas. 

Although it is possible that the doublet structure is due 
to the splitting of the exciton transition, as described by 
our model, there are certain problems with this inter- 
pretation. Firstly, it is not clear that one should even 
observe any measurable signal from the higher energy ex- 
citon states: these states lie more than 0.5meV above the 
lowest energy exciton transition, so, if the exciton were 
to be in thermal equilibrium, the expected populations 
of these states would be vanishingly small. Secondly, as- 
suming that there is some non-equilibrium population of 
the higher energy exciton states, one would perhaps ex- 
pect to see more than one additional line (recall that for 
a strongly-confined interstitial exciton all excited states 
have the same luminescence strength). Furthermore, this 
is not the only possible explanation for this structure. For 
instance, the recombination from a negatively charged 
exciton, which has been resolved in low-density GaAs 
wells at weak field p5|,^ and in CdTe quantum wells in 
strong magnetic field |27f|, may also be expected to appear 
in these systems, and has not as yet been identified. We 
are therefore cautious to ascribe the observed doublet to 
the line splitting predicted by our model, which may not 
appear in these experiments due to a lack of population of 
the higher energy exciton states. Further experiments are 
required in order to identify this structure. In particular, 
a study of the optical absorption spectrum would be very 
valuable in this respect. The exciton transitions that we 
have discussed will present much stronger spectral fea- 
tures in this experiment than in photoluminescence since 
the population factors, which reduce the observed inten- 
sities of high-energy exciton states in photoluminescence, 
are removed. 
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VI. SUMMARY 

We have presented a model ol photoluminescence in 
the presence of a magnetically-induced Wigner crystal, 
for systems in which the photo-excited hole lies close to 
the interface. This model predicts a splitting of the ex- 
citon transition due to the scattering of the exciton by 
the electrons forming the Wigner crystal. We studied the 
behaviour of our model as a function of the separation 
d. For d < I, the exciton is spread over many unit cells 
of the crystal. Additional peaks with small luminescence 
strength can appear in the photoluminescence spectrum, 
arising from higher energy states in which the exciton 
has a wavevector equal to a reciprocal lattice vector of 
the crystal. The spectrum carries diffraction information 
and is therefore sensitive to long-range crystalline order. 
For d > I, the exciton becomes strongly confined to an 
interstitial region. We studied this limit by developing 
an effective theory for the centre-of-mass motion of the 
exciton. Within the harmonic approximation for the po- 
tential at an interstitial site, the spectrum consists of an 
equally spaced set of peaks with relative intensities deter- 
mined by the relative populations of the various states. 

We discussed the photoluminescence experiments of 
Ref. [D and showed that the energy splitting of the dou- 
blet which has been associated with the presence of a 
Wigner crystal is comparable to the energy splitting we 
expect from our model for d ~ 100 A. We argued that this 
could be due to the recombination from a interstitial ex- 
citon. However, other explanations are also possible, and 
further experimental investigation is required to identify 
the structure our model predicts. In particular, this could 
best be observed in optical absorption spectroscopy. 

Whilst preparing this work for publication, we learned 
of work by J. R. Chapman, N. F. Johnson and 
V. N. Nicopoulos, proposing an interpretation of the pho- 
toluminescence experiments that is somewhat different 
from that presented here ||2^; we are grateful to them 
for a preprint of their work. 
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